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The Green-function formalism for the electromagnetic field in a 
magnetoelectric (ME) medium is constructed, as a generalization of 
conventional Casimir theory. Zero temperature is assumed. It is 
shown how the formalism predicts electromagnetic momentum to be 
extracted from the vacuum field, just analogous to how energy is ex- 
tracted in the Casimir case. The possibility of extracting momentum 
from vacuum was discussed recently by Feigel [Phys. Rev. Lett. 92, 
020404 (2004)]. By contrast to Feigel's approach, we assume that the 
ME coupling occurs naturally, rather than being produced by external 
strong fields. We also find the same effect qualitatively via another 
route, by considering one single electromagnetic mode. 
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1 Introduction 

Consider a magnetoelectric (ME) medium whose constitutive relations can 
be written on compact form as 
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Here x is the ME coupling parameter, assumed in general to be a pseu- 
dotensor, with = Xki- We employ SI units, so that the relation 

Sqixq = Ijc? refers to a vacuum, and we let the permittivity tensor e^^ and 
permeability tensor [i^^ be nondimensional (i.e., relative), so that the rela- 
tions Di = EoEikEk and Bi = iiQixikHk apply in the non-chiral case when 
X = 0. We shall take all material quantities Eik, fiik,Xik to be real and fre- 
quency independent. The tensors and fitk are symmetric; this being a 
general electrodynamic property following from the symmetry of the kinetic 
coefficients [Ij. No such symmetry condition exists for x, however. In some 
materials Xik is symmetric, Xik = Q^ik with g a pseudoscalar function, or, 
what is of more interest in the present context, Xik can be antisymmetric. An 
anisotropic crystal is called biaxial if the diagonal permittivity components 

^ £y ^ along the principal axes, and is called uniaxial if Sx = Sy ^ £z- 

In the following we will focus attention on the situation where the anisotropy 
in Xik occurs naturally. Cases where the anisotropy is created artificially, by 
means of strong electric and magnetic fields perpendicular to the direction 
of light propagation (cf., for instance, Ref. [2]), are for the most part outside 
the scope of the present paper. 

The macroscopic theory of ME media has been known for a long time. 
The reader may consult the book of O'Dell, for instance [3], as well as classic 
papers [HE]. A recent review is given by Fiebig [6] ; other relatively recent 
papers are Refs. [3 [HI El HO]- As explained in the Fiebig paper, two major 
sources for "large" ME effects can be identified: (i) In composite materials 
the ME effect is generated as a product property of a magnetostrictive and 
a piezoelectric compound. A linear ME polarization is induced by a weak 
ac magnetic field oscillating in the presence of a strong dc bias field, (ii) 
In multiferroics the internal magnetic and/or electric fields are enhanced by 
multiple long-range ordering. The ME effect can be strong enough to trigger 
magnetic or electrical phase transitions. 

The recent paper of Feigel [11] - cf. also the comments [12], [131 [HI [IS] - 
sharpened the interest in this special kind of materials. The main idea of this 
paper was to suggest a new quantum mechanical effect, namely the extraction 
of material momentum from the electromagnetic vacuum oscillations. The 
suggested effect is thus analogous to the well known Casimir effect [H], in 
which case it is an energy, not a momentum, that is extracted from the 
vacuum field. The Feigel effect thus belongs to a very active area in modern 
physics. Its main theme is the observability and the interpretation of vacuum- 
induced phenomena in macroscopic media. The effect has moreover a bearing 
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on the famous Abraham-Minkowski energy-momentum problem in dielectric 
matter [T7l[T5]. 

And this brings us to the main topic of the present paper, which is to 
investigate how the Green function approach, frequently used in Casimir- 
related problems, can be applied to a ME medium. To our knowledge, such 
a general approach has not been developed before. We follow the same basic 
field theoretical method as in the recent paper of Ellingsen and Brevik [12] , 
dealing with the Casimir effect. We will show that, even in the presence of 
the complexity in formalism caused by the ME effect, the theory leads to a 
right/left asymmetry in a medium- filled cavity enclosed within conducting 
walls placed at positions z = and z = a, and thus permits the extraction 
of momentum from the vacuum field, in principle. Our field theoretical for- 
malism thus supports earlier results that were based upon consideration of 
particular modes only. We will also have the opportunity to comment occa- 
sionally on some of the papers that followed the Feigel paper [201 HU |22] . 

In sections 2 and 3 we establish the governing equation for the Green 
functions, relate this to the two-point functions for the electromagnetic fields, 
and give explicit solutions in the presence of the two conducting plates. In 
Sect. 4 we digress to consider the momentum conservation equation for a 
ME medium, and show how the right/left momentum asymmetry occurs for 
one single mode. In Sect. 5 we return to the Green-function approach, and 
show how the momentum asymmetry occurs also for the vacuum field, when 
summing over all modes propagating in the ±x directions. 

We thus discuss the momentum asymmetry via two different approaches. 
A more detailed overview of the outline of the paper is given in Sect. 6. 

Readers interested in recent reviews on the Casimir effect may consult 
Refs. [231 El [23 ESI [23 EB]. Much information can also be found in the 
recent special issues of J. Phys. A [2D] and New. J. Phys. [5U] . 

We emphasize that the formalism below is constructed from the same 
main standpoint as in conventional Casimir theory: we calculate the change 
in field momentum caused by the geometric boundaries, i.e., the plates. The 
undisturbed system with respect to which we regularize Green-function ex- 
pressions is an infinite medium (without plates), made up of the same mate- 
rial. It is thus clear that in the limit when the separation between the plates 
goes to infinity, the effect that we calculate has to go to zero. 
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2 Governing equations for Green's function 



In this section we will establish the governing equations for the retarded 
Green function in the chiral medium. When this function is known, one can 
find the electromagnetic two-point functions and thus construct expressions 
for energy and momentum in the field. From now on, we assume the material 
to be isotropic, so that Eik = eSik, i^ik = l^^ik- Important in our context is 
that the coupling parameter Xik will still be permitted to be anisotropic. As 
already mentioned we take all material parameters e, /i, Xik to be real and 
frequency independent. They will moreover be assumed to be independent 
of the spatial coordinates. Our medium is thus assumed to be spatially 
homogeneous but chiral. (If the anisotropy of Xik is created artificially, by 
means of strong crossed electric and magnetic fields, the anisotropy property 
of Xik holds of course only in the constant field region between the condenser 
plates.) 

Let us first invert the constitutive relations ([T]) and ([2D to get 

E = — (d-^X-h), (3) 

ESQ \ C J 

B = /i/io(H + ^x'^-D). (4) 

These expressions hold when the EM effect is small, \xik\ ^ 1, what in 
practice always is the case. Terms of order x^ are neglected. 
Consider now Maxwell's equations in conventional form 

V ■ D = p, V ■ B = 0, (5) 

VxE = -B, VxH = J + D, (6) 

and take the curl of the first member of ([6]). Observing Eq. (j4]) we then get, 
when neglecting terms of order throughout, the following coupled vector 
equation for the basic fields E and B 

VxVxE + ^E + ^X-B + -Vx (x^-E) = -/ipoJ- (7) 

c c 

If X = 0, the coupling between the fields is absent. On component form the 
equation can be written 

V^-Ej - (9i (V ■ E) - - -XikBk - -Xik curlikEi = fifioJi. (8) 

c c 
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We have here defined curia. = eij^dj, where eijk is the antisymmetric Levi- 
Civita symbol with eijk = 1. 

In Eq. ([8]), the magnetic field Bk can actually be replaced by electric field 
components in view of one of Maxwell's equations, Bk = —curlkiEi. We 
obtain 



6uV^ - didi - 6ii^ dt + -Xik curlu dt - -Xik curkk dt 

c c 



El = fifioJi, (9) 



with dt = d/dt. 

We now turn to the Green-function approach. According to the source 
theory of Schwinger et al. (see, for instance, Refs. [25 or [31]), we make the 
correspondence J ^ P. p — —V ■ P. We introduce a dyad r(x,x') such 
that 

E(x) = — / c/Vr(x,x') ■ P(x'), (10) 

where x = (r,t). Due to causality, t' is only integrated over the region 
t' < t. The dyad F is the retarded Green function; also called the generahzed 
susceptibility. We take the Fourier transform of F, 

F(x, x') = r ^ e— F(r, r',co), T = t- t', (11) 

exploiting the stationarity of the system. We transform also the electric field, 

E(x)= r|^e-*E(r,c^), (12) 

with a similar expression for P(x). The governing equation for the Green 
function then becomes 

V X V X F(r, r', uj) - ^^F(r, r', lu) + — x ■ [V x F(r, r', lu)] 

& c 

- X [x^.F(r, r', .;)] = ^5{r - r')l, (13) 

c 

or, on component form. 



didj - dijV —dij H Xii curkj 



& c 

—XjiCurlu^Tjk{r,r,uj) = o{r - r )dik. (14) 



If Xik = and /i = 1, this equation reduces to Eq. (75.16) in Ref. [22] (their 
symbol is the same as our —hc^Tik/uj'^). 

If rik{r,r',u)) is known, we can make use of the fluctuation-dissipation 
theorem (which has a meaning both classically and quantum mechanically; 
cf. Refs. [321 [33]), to calculate the two-point functions: 

t{E,{r)Ekir'))^ = -Im{r,,(r, r', cu)}, (15) 

i{Bi{r)Bk{r'))^ = —\ curkjcurl^^ lm{Tji{r, r', u)}, (16) 
{Ei{r)Bk{r'))^ = --cwr^ Im{r,Kr, r', ^)}. (17) 

So UJ 

Here curl[j^ = eijkd'j, where dj is the derivative with respect to component j 
of r'. The expressions (fT5l) - (fT711 refer to zero temperature; a factor sgn(a;) is 
omitted throughout. The spectral correlation tensor (i?j(r)i?fc(r'))^ is defined 
according to 

— e~^-^{E,{v)E,{v'))^. (18) 

-oo ZTT 

[Note the meaning of the formalism here: the spectral correlation tensor is 
related to the Fourier transform {Ei{r, uj)Ek{r', uj')) of the two-point function 

{E,{x)Ekix')) via 

{E,{r,uj)Ek{r',uj')) = 2ti{E.,{v)Eu{v'))J{uo + uo')- (19) 

cf. Eq. (122.12) in Ref. |33j or also Appendix B in Ref. [3l].] 

Before going on to solve these equations, we will specify the geometry to 
be assumed in the rest of this paper. 

3 Specification of the geometry. Solutions for 
the Green functions 

Let us assume the same setup as in conventional Casimir theory, namely two 
perfectly conducting parallel plates separated by a gap a. The geometry is 
sketched in Fig. 1. 
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Figure 1: Sketch of the geometry. The TM and TE modes are shown. The 
tensor Xik^ constant everywhere in the fluid, is given by Eq. fl20l) . The wave 
vector k is directed along the x axis. 
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As mentioned earlier, we will mainly be considering the case where the 
ME effect occurs naturally. We assume accordingly that Xik is given initially 
and is constant everywhere in the ffuid, on the inside as well as on the outside 
of the plates. Because of the translational invariance in the x and y directions 
we can transform the Green function once more to obtain 

r(r,r',c.) = J -0^^e^^<^-^'^g{z,z',Kuj). (20) 
We also transform the delta function: 

5(r-r') = /^e^'^-(-')5(z-/), (21) 
and assume that Xik has the following form: 



Xik 



/ 
X,. I (22) 

\0 Xzy 



Our coordinate system is thus henceforth fixed, relative to the material. 
We focus attention on only one particular wave number k in the following, 
namely k = k^e^, directed along the x axis. In Fig. la) and b) the trans- 
verse magnetic (TM) and the transverse electric (TE) modes in the cavity 
corresponding to this k vector are indicated [55] . 

Our conventions above mean that we can let —>■ d"^ — k^. We can now 
write down the governing equations for the Fourier components Qik, from 
Eqs. ( IT1|) . The simplest equation follows by setting [ik) = (yy): 

9z9yy - (« — Xyz)9yy = ^S{z ~z), (23) 

where we have defined 

= kl - eiiuj^/c\ (24) 

Equation fl2^ is uncoupled; this being a consequence of our choice for k 
implying that dy ^ 

Setting (ik) = (xx) we obtain 

2 2 

i{kx + ^Xzy)d,g,x - {dl + ^^)fc = ^^{.z - z'), (25) 
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and with (ik) = (zz), 



i{kx + —Xzy)dz9xz + (/«^ + '^^^—^ Xzy)9zz = ^ 5{z - z'). (26) 

The last two equations are coupled. Consider finally the nondiagonal com- 
ponents: with (ik) = (zx) we obtain 

i{kx + ^Xzy)dzgxx + (/t^ + '^^^J^ Xzy)gzx = 0, (27) 

and with (ik) = (xz), 

^{kx + ^Xzy)d.9zz - {dj + '-^)9xz = 0. (28) 

The coupling in the differential equation (l25l) for Qxx can be removed if we 
make use of Eq. (!27|) differentiated with respect to x. Some manipulations, 
again observing that Xik is small, yield 



dlgxx-K'gxx = —6{z-z'), K = ^{l + ^Xzy). (29) 



Equation fl23|) can be rewritten similarly: 



9yy - L'^9yy = ^^K^ " ^')> L = k{1 - ^^^^Xyz)- (30) 



The differential equations ( 1291) and (1301) for the diagonal components are 
convenient for further manipulation. Note that the values of K and L are 
dependent on whether the direction of propagation of the wave is to the right 
or to the left. If Xik = 0, the expressions agree with those of Ref. [19] [36] . 

We now proceed to solve the equations, beginning with Eq. ( l30l) . As 
Ey = aX z = and z = a because of the boundary conditions, we have 
gyy{0, z' ,'k,u) = gyy{a, z' ,k,u!) = 0. The solution of Eq. ( l30l) can then be 
written 

_ j^^L\z-z'\ __ ^~L(z+z') , 2 [cosh L{z - z') - cosh L{z + z')] \ 
2Lc^\ exp(2La)-l j' 

(31) 

When Xik = 0, this expression agrees with that given in Appendix C of 
Ref. [37] in the limit of perfectly conducting plates. 
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Next considering g^x, we must analogously have QxxiS^i z' , k, oj) = gxx{ci, z', k, 
in view of the boundary conditions. The solution of Eq. (1291) becomes 

[ j^l,„,/| r^(^^^r) 2[cosh K{z- z') -cosh K{z + z')] ] 

24 + exp(2Ka)-l j' 

(32) 

The expressions f l3T|) and fl32l) are fairly complicated. For practical pur- 
poses it is possible to simplify the expressions considerably, by omitting terms 
containing [z + z'). The reason is that these terms do not contribute to phys- 
ical quantities like the Casimir force on the plates or to the field momentum 
in the gap. This can be seen in two different ways. The simplest way is to 
argue, as in Sect. 81 in [32], that by putting z = z' in solutions having the 
argument {z + z') one would obtain physical quantities like field momentum 
in the gap varying with the position z. This would contradict the law of con- 
servation of momentum. Another way of examining this rather subtle point 
is to include the {z + z') terms everywhere in the formalism, and to verify 
that they really do not contribute in the end. In addition to the discussion 
in |32j, one can find more mathematical details about this point in the paper 
[in] and in the thesis [38j. 

We can moreover omit the source-dependent inhomogeneous \z — z'\ term 
in each of the Green functions. This term represents the solution pertaining 
to the delta function source inside a homogeneous medium filling all space. 
Being geometry independent, it cannot contribute to any physical quantity 
related to the geometry. All in all, we shall in the following use the "effective" 
Green functions 

_ K cosh. K{z - z') 
e exp (2a a) — 1 

_ /io^ coshL(2: - z') 

Lc2 exp(2La)-l- ^ ' 

Consider finally the remaining diagonal component, g^z- To this end we first 
observe the symmetry property 

gxz{z, z\ k, u) = gzx{z', z, -k, uj), (35) 

which is an example of the general relation 

Tikir, r', r) = Tki{r', r, -r), (36) 
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when expressed in Fourier space (cf. Sect. 81 in Ref. [52] )• From Eq. (1271) 
we have, when inserting the expression (155]) . 



/ /I \ "^/i /^^ , sinh (2; — 2;') . 
^..(^, . , k, = -(A:. + -X.,) ^^p(2ira)-l - 

Now, to the required order K^—k^) = K^/K{kx), according to Eq. (p9|) . Thus 
we get from Eq. (153]) 

i /iu; sinh[/t2(z - 

£ c exp(2K^a/A) — 1 



where here and henceforth K = K(kx) as defined in Eq. (l29l) . From Eq. (l26i) 
we then finally get (delta-function omitted): 

^ ^ cosh[/.^(^-^0/ir] 

7^3^ exp(2K2a/ir) - 1 ' ^ > 

These Green-function expressions have to our knowledge not been derived 
before. Before applying them to the Feigel effect, we shall in the next section 
follow a more simplistic approach and consider the right/left asymmetry in 
the field momentum considering one single mode only. 



4 Energy- momentum formalism. Right/left 
field momentum asymmetry 

4.1 Energy- momentum formalism 

Before considering the momentum asymmetry for one single chosen direction 
of propagation, we need to develop the formalism related to the electro- 
magnetic energy-momentum tensor. In this subsection we take a general 
approach, allowing for external charges p and currents J. The coupling 
tensor Xik is allowed to be general (not necessarily of the form given in 
Eq. (I22I) ). though constant, and we allow also for optical anisotropy by let- 
ting e5ik — > Sik, pSik — >■ Pik with the material parameters constant. 

It is convenient to write the constitutive relations ([1]) and ([2]) on tensor 
form, 

A = SoSikEk H XikBk, (40) 

/ioc 
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Hi = — —Xki Ek + —fi \k Bk- (41) 

/ioC /io 

In view of Maxwell's equations ([5]) and ([6]) we obtain the conservation equa- 
tion for energy, 

V-S + w = -E-J, (42) 
where E ■ J is the energy dissipation, 

S = E X H (43) 

the Poynting vector, and 

«; = -(E-D + H-B) (44) 
2 

the energy density. 

As for the momentum conservation, it is convenient to start from the 
equation 

<9t(D X B)i = —pEi — EijkJjBk + dk{EiDk + HiBk) — SkiEk^iEi — /i ^^iBi^iBk, 

(45) 

which follows from Maxwell's equations (here Ei^k = dkEi, etc.). Introducing 
the Lorentz force density 

= pE + J X B, (46) 
as well as the Minkowski stress tensor |17j . 

= E,Dk + HA - ^(E ■ D + H • B), (47) 
we can write the momentum conservation equation as 

- r = , (48) 

where 

g^^ = D X B (49) 

is the Minkowski momentum density. (The symbol g for momentum is not to 
be confused with the Green functions.) It is generally known that the above 
expressions hold when the medium is optically anisotropic. It is however 
somewhat remarkable that they hold when Xik 7^ also; there seems to be 
no simple physical reason why Xik should drop out from the formalism. 
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In the case of high frequency fields, in particular optical fields, the Minkowski 
theory appears to be both simple and capable of describing all experiments 
(cf. the analysis of one of the present authors on this point some years 
ago [39]; some more recent papers are listed in Ref. llSj). However, at low 
frequencies where the effect of the oscillations are themselves observable - 
notably in the Lahoz- Walker experiment [10] - the experiments agree not 
with the Minkowski but rather with the Abraham force, which accordingly 
can be taken to be the most 'physical' alternative at these frequencies. The 
Abraham theory [18] consists in symmetrizing the stress tensor, 

Tf, = ^{E,Dk + EkDi) + ]^{HiBk + HuBi) - ^(E ■ D + H ■ B), (50) 

and taking the momentum density to be 

= ^E X H, (51) 

the latter satisfying the relation g = S/c^, the so-called Planck's principle of 
inertia of energy. 

We assume henceforth optical anisotropy so that e and yU are scalars, 
and also that p = 0, J = 0. The Minkowski and Abraham stress tensors 
become thereby equal, = T^. The Abraham conservation equation for 
momentum can be written as 

dkT^^ - gt = [{e^^ - l)/c']dt{^ x H)„ (52) 

where the term on the right hand side is the 'Abraham term'. It was precisely 
this term that was measured by Walker and Lahoz [IQ]. In a high-frequency 
field, it fluctuates out. 

We shall return to the Abraham force in Sect. 6. 

4.2 Momentum asymmetry 

Referring to Fig. 1, we consider to begin with only the right-moving TE wave 
corresponding to the field components 

Ey = ^ sin knz e*(k-x-c.i) ^ ^53^ 
i?. = ./^^cosA:„^e^('^— (54) 
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5, = W-^sinfc^ze^C^-"-"*) 
y a u 



(55) 



the other components being zero. (We use the same normahzation of the 
fields as Tiggelen et a/. [22].) Here k-x = kxX, and kn = nn/a with n = 
1, 2, 3... is the transverse wave number. For a given value of k^, the eigenfre- 
quencies uj are thus discrete. We can derive the dispersion equation by going 
back to the field equation ([9]) for Ei = Ey in the source-free case, observing 
that diEi = dyEy = 0, inserting the form (122!) for Xik- We obtain 



kl + kl-'J^co^-^J^ 
\ c 

which implies to the lowest order in Xyz 



Xyz ) Ey 0, 



(56) 



1.2 I 1.2 



1 



^ x kl + kl 



Xyz 



(57) 



The right/left asymmetry is manifest. A left-moving wave is described by 
the substitution k^ —k^. 

Let us now calculate the field energy density, w, for the TE mode. We 

get 



w 



fE D* + H B* 



£o£ 
2a 



1 + 



SIV? knZ + 



kl 



COS IxijiZ ^ 



(58) 



2fiQfia up' 

the Xyz terms drop out when w is written in this way. It is convenient to 
consider the expression integrated from 2; = to 2; = a, thereby getting the 
energy W per unit length and width. 



W 



waz = — — 
4 



(59) 



Using Eq. (|57|) we can write this in terms of the wave number components. 



1 + 1 //^ 

2V £ Jk-^ + kl 



Xyz 



(60) 
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The Poynting vector in the x direction, Sx-, may be calculated as 



Sx = -(E X H* 



Xyz 



sin knZ, 



(61) 



which means that the integrated energy flux when expressed in terms of wave 
number becomes 



r c J IT kx 
qx= Sxdz = — J j= 



+ kl kl + kl 



Xyz 



(62) 



Alternatively, we might calculate the energy flux as qx = Wux, where Ux is 
the group velocity 



Ux 



duj 
dkr 



kfr 



Xyz- 



(63) 



This agreement is as we should expect, since we are dealing with the prop- 
agation of low-amplitude waves. The kinematic group velocity concept and 
the dynamic energy flow velocity concept should be the same. 
Consider finally the Minkowski momentum density g^: 



9^ 



1 



SqE k, 



-(D X B*), = ^— [l + ^xyz] sin' knZ 



a u 



(64) 



Comparison between Eqs. ( l6TI) and ( l64l) shows that the relationship = 
{e^/(?)Sx-, known from conventional optics, does not hold when Xyz is differ- 
ent from zero. We also give the expression ( 16^ when integrated over z: 



SqE kx 



Jo 

Again, the right/left asymmetry is manifest 



o \ 1 + Xyz 

2 uj \ ecu , 



(65) 



5 The Green-function approach to the Feigel 
effect 

Our intention now is to calculate the Minkowski momentum asymmetry in 
the chiral medium using the Green-function approach from Section 3. We 
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start from the following general expression, reverting to real representation 
for the fields, 



g 



M 



lim 



r V^'''^^ I ^2 ^""^'-'^ (D(r) X B{r')U (66) 



(27r 



We assume zero temperature, so that the brackets ( ) mean purely quantum 
mechanical average. As no thermal fields are excited, the field momentum as 
well as the field energy stem exclusively from the vacuum zero-point oscilla- 
tions. Whereas in the previous section we considered the contribution from 
one single selected mode only, we shall now consider the effect of summing 
over all available vacuum modes. We shall impose one restriction, however: 
the wave number k will be required to lie either in the positive or the nega- 
tive X direction. This corresponds to our Green-function approach in Sect. 3. 
Mathematically, it means that we can let / (P'k/{2TxY J dkx/27c. As the 
distribution of fields does not vary in the transverse y direction, we can ef- 
fectively let dy ^ when applied to the fields. Evidently, the x component 
of field momentum has to be zero in the case of a non-chiral medium; if there 
is an asymmetry following from the formalism this has to be caused by the 
presence of Xik- As before, we assume the particular form (!22l) for Xik- 
The X component of Eq. fl66|) becomes (we omit the 'lim' from now on) 

= II II '^^^''^^^^^'\Dy{v)B,{v') - D,{v)By{v')U. (67) 



We insert from Eqs. ([T]) and ([2]) 



Xyz 



Dy = EQeEy H Bz, 

/ioc 

Dz = SoeEz + ^^By, 
Hoc 



(68) 
(69) 



and get 



9^ 



271 



2n 



eoe{Ey{r)Bzir'))^k 



- eoe{Ez{r)By{r')U - ^{By{r)By{r'))^, + ^(i?,(r)5,(r')).fc . (70) 



Xyz 
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We have thus so far expressed in terms of the two-point functions for the 
fundamental fields. Using Eqs. (ITB]) and (fTTI) we calculate 



Eq uj 

{E^{r)By{r'))^k = — -(<9^ ^t^9zx + iK l^Qzz), 

60 to 

h i 

{By{v)By{r'))^k = ^(<9^ Im^ia;^ - ika: Imd^g^x 

- ikxdz Im - kl Im g^^), 
(5^(r)5^(r'))^fc = T^^T^9yy 

We have here, as above, naturally defined ( )u}k via the relation 



etc. We can thus express g^ as 



9x 



*^ n 



27TUJ 



/-oo 27r 



{x-x') 



where 



( ) = -ikxE Im gyy - E{d[ Im g^^ + ikx Im g^z) 
zc 

H — Im gxx + ^A^z Im d^g^^ 
to 

tk^ c 

+ ikxdz Im gxz + fc^ Im g^z) —Xyz Im gyy. 

LiJ 



(71) 
(72) 

(73) 
(74) 



/ 00 rjU 
-Le^M^-^){Ey{v)Bz{v'))^,, (75) 
~oo ZTX 



(76) 



(77) 



This expression shows that it is necessary to calculate gyy,gzx and gzz to 
order Xik- From Eqs. flMl) . fl37|) and fl39|) we get 



9yy 

dz9zx dz9zx 



nd 


1 + 


fJ'kxOO 




1 + 




Ed 1 


kxC 



1 + 



2/ta 



,2Ka 



! + %(! 



2Ka 



(78) 
(79) 
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9z 



end 



,2rea ' 



2Ka e 



Xzy 



where d is defined as 



d = e'^"" - 1. 



(80) 



(81) 



The remaining terms in Eq. (|77ll are however multiplying Xzy or Xyz, and so 
need not to be expanded in Xik- Thus to sufficient accuracy 



dzQzx dzQxz 



1 



1 

d' 



9yy 



9zz = ;• 



(82) 

(83) 
(84) 
(85) 



KC^ d 

Ml. 

ne d 

We now put t = 0, x — x' = in Eq. fl76|) . and perform a standard complex 
frequency rotation whereby u —>■ i(, with ( real ^31j. As du/co —>■ d(/(, it 
follows from Eq. fl751) that of physical importance are only those terms in ( ) 
that are real after the rotation {g^ has to be real). Thus the ffist terms in 
Eqs. ( !78|) . (1791) and (IHOj) do not contribute. This is what we should expect: the 
asymmetry in momentum is caused by Xik- After some calculation we obtain, 
by letting /f^ d( ^ 2 /g°° d(, dk^ 2 /q°° dk^ because of symmetry of 
the integrand about the origin, 



9^ 



Ahfj d( 


poo 


c Jo 


277 


Jo 


P~ 


1 + 


3kl 
2k' 



271 K^d 




klc' 



Xyz 



Xzy 



(86) 



Recall that d is given by Eq. (IHTj) . where now k' = k' + e[iC^I(?. The 
integrals are seen to be finite. This is so because we have already performed 
the regularization by omitting those parts in the Green function that refer to 
the infinite undisturbed system. (Cf. also the remarks at the end of Sect. 1.) 



If the separation becomes infinite, then d — ^ e 



2Ka 



OO. 



and 



M 



as it 



must; all plate- induced physical effects have to go away in this limit. 
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The expression fl86p may be conveniently rewritten in terms of polar co- 
ordinates. Introduce X = k^: = kcosO, Y = {^feJi/c)C, = nsinO, so that 



(87) 



The area element in the XY plane is ndndd = {y/6jl/c)dkxd(. Then 



,f = -^-y^ cosHdel ^|[l-tan^^— e-Jx.. 



-5 + 2 tan' 6 - {2 + tan' 6) Xzy 



(88) 



The integrals can be evaluated to give 




has) fji 



Xzy, 



(89) 




where ^(3) is the Riemann zeta function with argument 3. It is noteworthy 
that only one of the ME coefficients, Xzy, appears in this expression. The 
factor multiplying Xyz in Eq- (IHH!) turns out to be zero. There seems to be 
no simple reason for this, although the behavior is obviously related to the 
complicated structure of Eq. (170|) and the need to expand Qyy^Qzx and g^z 
in Xik] cf. Eqs. (rr5]) - (|5U]l . Recall that the x direction has been singled out 
as special, and also that we have taken the variation of the fields in the 
transverse y direction to be equal to zero. To our knowledge, an expression 
of this kind has not been derived before. 

The quantity g^^ is measurable, in principle. Before any measurement can 
be done, the expression ( l89l) has of course to be augmented by contributions 
from all the other values of k. Aspects connected with real experiments lie 
outside the scope of the present paper. 

6 Summary, and discussion 

Let us first recall the assumption on which the above calculation is based: 

1. The tensor Xik characterizing the magnetoelectric medium is given 
naturally over all space in the fiuid, on the inside as well as on the outside 
of the conducting plates. The constitutive relations are Eqs. ([1]) and ([2]), or, 
in inverse form, Eqs. ([2]) and (jlj) when the magnitude \Xik\ of the coupling is 
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small. In the example that we calculated in detail, Xik is given by Eq. fl22p . 
The tensor Xik may be asymmetric, in contrast to the permittivity Sik and 
permeability fiik which are always symmetric. 

2. We have followed two different approaches, giving most weight to 
the Green-function approach since this does not seem to be treated very 
much in the literature. We took the temperature to be zero. The governing 
equation for the dyad Fj^ is Eq. f lTll) . The full solution of two of the diagonal 
components, Qyy and Qxx, introduced as Fourier components of the F's via 
Eq. (l20l) . are given by Eqs. (!3T|) and (!32|) . We have here assumed that there 
is no variation of the fields in the transverse y direction. The derivation of 
Eqs. (13 ip and (1320 generalizes conventional Green-function Casimir theory 
[2H [57] to the case of ME media. For practical purposes it turns out to 
be possible to simplifying the expressions considerably, by omitting terms 
that do not contribute to physical quantities in the end. The arguments for 
proceeding in this way are spelled out, for instance, in Ref. [32]. The relevant 
reduced components of Qik in our case are given at the end of Sect. 3. 

3. In Sect. 4 we deviated to follow a different, and more simple, approach. 
After having established the momentum conservation equation for a ME 
medium, we calculated the right/left asymmetry for one single mode only 
(choosing one of the modes considered in Ref. p2]). The results are given by 
Eqs. (!MI) and (l65l) . Adding two similar modes, one propagating in the +x 
direction and one in the —x direction, we obtain a net flow of momentum, 
caused by the coupling Xyz- 

4. In Sect. 5 we returned to the Green-function approach, calculating the 
net X component of momentum arising now not from one single mode, but 
from all modes propagating in the ±x directions in the vacuum field. The 
main result is given by Eq. (15^ . All terms independent of Xik drop auto- 
matically out of the formalism, in accordance with what we should expect 
beforehand. 

5. On physical grounds one may ask: where does the net electromagnetic 
momentum come from? Obviously, it cannot come from 'nothing'. We are 
actually comparing two different physical situations here. The first is when 
the conducting plates are infinitely far separated. This is our initial 'vacuum' 
state. The final state is when the plates have been brought close to each other, 
infinitely slowly. The calculated quantity g!^^ is the Minkowski momentum 
density extracted during this process of change of the plate separation. The 
coupling parameter Xik in the fluid is the same, all the time. The process 
is thus conceptually quite close to the process encountered in usual Casimir 
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theory; the main difference being that it is now momentum, not energy, that 
is extracted. 

6. The setting of our thought experiment is similar, but not exactly the 
same, as that envisaged in Feigel's paper yjj. Feigel assumed the coupling 
Xik in the fluid to be the result of applying strong electric and magnetic fields. 
We have deliberately avoided this picture since it complicates the situation 
in the sense that one has to deal with two sets of fields, both the external 
fields, and the wave modes. When assuming naturally occurring Xik instead, 
as we have done, the interpretation of the effect becomes more transparent. 

Before leaving this idea, let us however not the following point: Assume 
that strong crossed fields Eq and Hq are applied between the conducting 
plates at the instant t = 0. Then, during the time when the external fields 
increase in strength, there acts an Abraham force in the fluid in the interior. 
The force density is given by the expression on the right in Eq. fl52|) . Inte- 
grating over time, from t = until the external fields have become constant, 
we see that the following mechanical momentum density is imparted to the 
fluid: 

g^ = ^^(EoxHo). (90) 

This is the dominant momentum given to the fluid between the plates. In 
addition comes the momentum transferred from the wave modes; these are 
connected with Xik- The momentum (189|) is actually very similar to the 
momentum (or more strictly the angular momentum) transferred to the sus- 
pended dielectric cylindrical shell in the Walker-Lahoz experiment in ordinary 
electrodynamics [40l [39]. 

7. It might appear surprising that in Feigel's paper a high-frequency 
cutoff ujcut is introduced, whereas in the present treatment there is no need 
of a cutoff. The reason for this behavior is that the two formalisms are 
constructed differently: Feigel considers the total contribution, including that 
of the infinite unconstrained system, whereas in our case we have regularized 
the infinite contribution away. The case of high frequencies leads in Feigel's 
case to infinities, whereas in our case it leads to zero. Again, this is the same 
point as was emphasized at the end of Sect. 1. We are generally looking at 
the present problem as a sort of Casimir-type problem. 

8. What is the connection between the Feigel effect and relativity? In 
this context it might be of interest to recall how the relativistic formulation 
of electrodynamics in continuous media is formulated. There is always one 
particular inertial system here, namely the one where the medium is at rest 
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- this was emphasized aheady in the classic papers of Jauch and Watson [H] . 
The relativistic formulation is obtained by introducing two electromagnetic 
field tensors F^^j and H^j^y such that the covariant Maxwell equations 

dpF^, + d^F.p + d,Fp^ = 0, d,H^, = (91) 

agree with the standard Maxwell equations in (we assume no external 
charges or currents). The electromagnetic energy-momentum tensor Sfj,^,, 
assuming Minkowski's expression for the momentum density, is divergence- 
free, 

d.S^, = 0, (92) 

meaning that the energy and momentum of the total field constitute a four- 
vector. Moreover, this four-vector is space-like, so that it is possible to find 
inertial systems where the radiation energy becomes negative. A striking 
demonstration of this property is found in connection with the Cherenkov 
effect, in the frame where the emitting particle is at rest. A clear introduc- 
tion to this kind of theory is found in M0ller's book [42j, and the theory is 
discussed also in papers of one of the present authors [321 113] • 

In our opinion there is no strong connection between the Feigel effect and 
relativity. The force on the fluid, or the momentum transferred to it, are 
calculated assuming the fluid to be at rest. Relativity is as little involved 
here as it is involved in the description of the Walker-Lahoz experiment. An 
exceptional case is, however, if the Euler-Heisenberg Lagrangian is drawn 
into consideration as a model to describe the ME effect (cf. for instance, van 
Tiggelen et al. [22]). 

9. It is of interest to have an idea about the magnitude of the effect 
that we have considered. Magnetoelectric birefringence is actually found 
even in a vacuum, when there are strong crossed external fields Eq and Hq 
present. The effect is however extremely small. Let An = hb — ue denote the 
difference in the refractive index between the magnetic and electric directions. 
Even with a strong magnetic field of 30 T and an electric field of 10* V/m 
the birefringence is only An ^ 8 x 10^^'^ [7J. 

A more promising case is when one applies strong orthogonal fields to a 
linear isotropic liquid. Thus Roth and Rikken [2] performed an experiment in 
which molecular liquids were placed in such a strong field region. By passing 
laser light through the liquid, perpendicular to the fields, they obtained a 
linear relationship between the field strength and the MR birefringence. With 
a magnetic field strength up to 17 T and an electric field of 2.5x10^ V/m 
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the ME birefringence was found to be of order An ~ 10~^^. Thus the ME 
effect is much larger in a hquid than in a vacuum. 

Naturally occurring anisotropies, the case that we have been considering, 
seem actually to be stronger. Thus the crystal FeGaOs is known to be mag- 
netoelectrically active with ME coefficients about 3 x 10~^ at low frequencies. 
In this crystal, as well as in analogous crystals like FeAlOs, anisotropies of 
order 10~^ are expected over a wide frequency range from DC to X-rays [H] . 
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